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A linearly polarized light normally incident on a semiconductor quantum well with spin-orbit coupling may 
generate pure spin current via direct interband optical transition. An electric photocurrent can be extracted from 
the pure spin current when an in-plane magnetic field is applied, which has been recently observed in the In- 
GaAs/InAlAs quantum well [Dai et al. Phys. Rev. Lett. 104, 246601 (2010)]. Here we present a theoretical 
study of this magnetoelectric photocurrent effect associated with the interband transition. By employing the 
density matrix formalism, we show that the photoexcited carrier density has an anisotropic distribution in k 
space, strongly dependent on the orientation of the electron wave vector and the polarization of the light. This 
anisotropy provides an intuitive picture of the observed dependences of the photocurrent on the magnetic field 
and the polarization of the light. We also show that the ratio of the pure spin photocurrent to the magnetoelec- 
tric photocurrent is approximately equal to the ratio of the kinetic energy to the Zeeman energy, which enables 
us to estimate the magnitude of the pure spin photocurrent. The photocurrent density calculated with the help 
of an anisotropic Rashba model and the Kohn-Luttinger model can produce all three terms in the fitting for- 
mula for measured current, with comparable order of magnitude, but discrepancies are still present and further 
investigation is needed. 

PACS numbers: 42.65.-k, 72.25.Dc, 72.40.+W, 73.63.Hs 
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I. INTRODUCTION 

Optical injection and detection of spin current have at- 
tracted a lot of interest recently in the context of non- 
magnetic semiconductor spintronics.'-^ Exciting progress has 
been made in recent years, such as the injection of spin- 
polarized currents by using the spin galvanic effect^ and the 
circular photogalvanic effect,'* the observation of spin accu- 
mulations induced by the spin Hall effect with the help of the 
Kerr rotation^ and the p-n junction light-emitting diode,'' the 
inverse spin Hall effect,^-^ and the detection of pure spin pho- 
tocurrent by the second-harmonic effect.^'" 

Generation and detection of the pure spin current where 
there is no net charge cuiTent is one of the challenging issues 
in spintronics. Optically, the generation of pure spin photocur- 
rent has been realized by injecting orthogonally polarized one- 
and two-photon linear lights," '^ or by the incidence of un- 
polarized or Unearly polarized lights into bulk III-V semicon- 
ductors and quantum wells with spin-orbit couplings.i^^ In 
the latter case, a linearly polarized light can be decomposed 
into circularly polarized lights with opposite helicities. Due 
to the conservation of angular momentum, the circularly po- 
larized lights with opposite helicities excite caiTiers with op- 
posite spins, which are locked to opposite momenta owing to 
the spin-orbit coupling. As a result, the photo-injected carri- 
ers with opposite spins always have opposite velocities, lead- 
ing to the pure spin photocurrent. One way to observed this 
underlying pure spin photocurrent is to apply an in-plane mag- 
netic field, which will create an imbalance of spins. The spin 
imbalance leads to the asymmetry of electron velocities by 
the spin-orbit coupling and can extract an electric photocur- 
rent from the pure spin photocurrent. The field-induced con- 
version from pure spin photocurrent to electric photocurrent 



(also referred to as the magneto-gyrotropic photogalvanic ef- 
fect, MPGE) was systematically studied via the intra-band22i21 
and inter-subband21 electron heating by THz and microwave 
radiations, in which spin-dependent scatterings are necessary 
for the observed charge current. 

Very recently, the magnetoelectric photocurrent has been 
observed via direct interband optical transitions.^^ Unlike the 
intra-band transition, in the direct interband transition the 
electron momentum is conserved, and the magnetoelectric 
photocurrent is a few orders of magnitude larger than that via 
the intraband transition. Therefore, it appears to be promising 
as a simple and practical method to generate and detect the 
spin current. Although symmetry analysis^^ can justify the 
observed electric photocurrent as a function of magnetic fields 
and polarization of the light, a microscopic description is still 
needed. In particular, the current induced by a parallel field is 
against intuition, considering the form of the spin-orbit cou- 
pling derived from the symmetry argument. Moreover, how to 
estimate the zero-field pure spin photocuiTent from the mag- 
netoelectric photocurrent is an interesting issue. 

In this work, we study theoretically the generation of pho- 
tocurrent via direct interband transitions excited by shedding 
a linearly polarized light normally into a semiconductor quan- 
tum well. We find that the k-space anisotropy of the carrier 
density provides an intuitive microscopic picture for the mag- 
netoelectric photocurrent and the zero-field pure spin pho- 
tocurrent. The origin of the observed dependences of the 
photocurrent on the magnetic field and the polarization of the 
light can be illustrated within this picture. The photocurrent 
density calculated with the help of a minimal model with an 
anisotropic Rashba model and the Kohn-Luttinger model is 
comparable with the experiment. Moreover, we propose an 
approach to estimate the magnitude of the undetectable pure 



spin photocurrent, which is generated at zero magnetic field 
by the same linear light, from the observed magnetoelectric 
photocurrent. Part of the results was very briefly reported in 
Ref.|25t 

The paper is organized as follows. In Sec. |ll] we review 
key features of the experimentally observed magnetoelectric 
photocurrent as a function of magnetic fields and the polariza- 
tion of the linear light, and present a symmetry argument. In 
Sec. Hn] a minimal model of the quantum well is presented, 
as well as the electric -radiation interaction that accounts for 
the interband optical transitions. In Sec. IIVI we introduce the 
anisotropic photoexcited carrier density in k space, with the 
help of a standard density matrix formalism. In Sec. IV Al we 
illustrate the origin of the field and polarization dependences 
of the magnetoelectric photocurrent in terms of the anisotropic 
photoexcited carrier density. In Sec. IV Bl we compare the 
calculated magnetoelectric photocurrent densities with the ex- 
periment. In Sec. [yll the estimate of the underlying zero-field 
pure spin photocurrent is discussed. Finally, a summary is 
given in Sec. I VIII 



II. EXPERIMENT AND SYMMETRY ANALYSIS 

In this section, we will review the experiment on the mag- 
netoelectric photocurrent generated via direct interband opti- 
cal transitions by shedding a linearly polarized light normally 
into an InGaAs/InAlAs quantum well.^^ The experiment setup 
reported in Ref. |25|is schematically illustrated in Fig. [T] 
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FIG. 1 : Schematic illustration of the setup. A linearly polarized light 
(downward arrow) is normally incident into a semiconductor quan- 
tum well grown along the [001] direction of zinc-blende materials. G 
denotes the light polarization angle with respect to the x-axis. The 
X, y, and z axes are defined along [lIO], [110], and [001] crystallo- 
graphic directions, respectively. 6, {By) indicates the parallel (per- 
pendicular) magnetic field with respect to the generated photocurrent 
A. 

In the experiment, the current density measured along the x 
direction can be formulated as^ 



jx - coBy + CyBy COS 26 + CyBx sin 29, 



(1) 



where B^ and B, are the magnetic fields along the x and y 
directions, respectively. 9 is the polarization angle of the lin- 
early polarized light with respect to the x axis, co/.v/y are con- 
stant coefficients linearly scaled with the light power The x 
and y axes here are defined along [110] and [110] crystallo- 
graphic directions of the zinc blende structure, respectively. 



The experimental results are consistent with the symmetry 
argument of the Civ point group. ^^ Inversion-asymmetiic zinc 
blende heterostructures grown along [001] crystallographic 
direction usually have the Civ point group symmetry. As 
a result, the electric photocurrent induced by in-plane mag- 
netic fields can be phenomenologically written as (refer to 
AppendixlAlfor details) 



Ji' 



Y^ix^'^'B^EyE- +x'''^"BpEyEy), 



(2) 



where ju, v run over (x, y), and yu, v = y if yu, v = jc, and vice 
versa. ^ is a fourth-rank pseudo tensor that relates the elec- 
tric photocurrent to the polarization electric field components 
{Ex, Ey) oc (cos 9, sin 9) of the incident light and the applied 
magnetic fields {By, By). One can readily check that Eq. (|2| 
yields the same form as Eq. ([1]). 




FIG. 2: (Color online) Schematic description of direct interband opti- 
cal transitions from a valence band |v, k) to the spin-split conduction 
bands \c = ±,k) at zero magnetic field (a) and in the presence of a 
magnetic field B (b). The states involved in the transitions form a set 
of constant energy contours (denoted by dashed rings), which can be 
shifted along ky {k^) in the presence of a magnetic field along x {y) 
direction. 

Another consequence of the C2,. symmetry is that when x\\ 
[110] and yll [110], the invariants of the Civ symmetry allow 
the spin-orbit coupling in linear k to be only in the form.^^i21 
(refer to AppendixlAlfor details) 



^SOC = ^vO-xky - AyO-yky 



(3) 



where ky^y is the wave vector along the x {y) axis, the Pauli 
matrix cr yjy depicts the spin along the x {y) direction, and /Iv/y 
is the spin-orbit coupling coefficient along the x {y) direction. 
Usually, the current along x direction coiTesponds to the av- 
erage shift of momenta of carriers along the x direction. Ac- 
cording to Eq. ([3j, ky couples only to the spin along y direc- 
tions, thus it can only be shifted by a magnetic field along the 
y direction. A magnetic field applied along x direction can not 
shift the average k^ [see Fig. |2|b)]. According to this argu- 
ment, the last term of Eq. ([T) seems against intuition because 
it originates from a magnetic field along the x direction. 

The main task of this work is to present an intuitive picture 
to solve the above dilemma, as well as the magnetic field and 
polarization angle dependences of the electric photocuiTent, 
in terms of the anisotropy of the photoexcited caiTier density 
in k space. We will see that the three terms in Eq. ^ can be 
well interpreted by the microscopic picture presented in Fig. 
|4](subfigures in boxes). 



III. MODEL HAMILTONIAN 

In this work, we consider a minimal model that describes 
the energy bands near the band gap A of the quantum well. 
The advantage of the minimal model is that it is analytically 
solvable, thus physically transparent in delivering the physical 
picture. We expect that more sophisticated calculations, e.g. 
by using the 14-band k ■ p model'^ or the full band structure 
local density approximation^^ can cover more details, but will 
not deviate from the physical picture obtained by the minimal 
model. 

The total Hamiltonian is given by 



H = Hc + Hv + V(t). 



(4) 



where He is for the conduction bands, Hv for the valence 
bands, and V{t) for the electric-radiation that couples the con- 
duction and valance bands by the linearly polarized light. 
The lowest conduction subbands are described by a two- 
dimensional free electron gas with the Rashba and Dressel- 
haus spin-orbit couplings and in-plane magnetic fields, 



He = ^(kl + kl) + (A,cr,k, ■ 



'ixO-ykx) + h ■ O", 



(5) 



where m* is the eff'ective mass of electron, fi is Planck's con- 
stant over 2n, k^ (ky) is the wave vector along the x {y) axis. 
Ay = a+/3, /If - a-/3, and a and/3 are the Rashba and Dressel- 
haus spin-orbit coupling coefficients, respectively. Through- 
out the work, x and y axes are defined along [110] and [110] 
crystallographic directions, respectively, cr - (cr^, o-y) is the 
vector of the Pauli matrices, h - (hx, hy) - ^geHsiBjc, Bv) is 
the Zeeman energies induced by the in-plane magnetic fields 
(B^, By), with ge the Lande ^-factor, and fiB the Bohr magne- 
ton. The eigenenergies and eigenstates of He are given by 



e„±(k) = -^k^ ± A, 
Zm 



l+,k) = 



1/V2 
U 



|-,k) = 



-U' 
1/V2 



(6) 



(7) 



where U - [Ayky + h^ - ((A^kx - hy)] / V2A and A - 

^{Ayky + /Z,)2 + {Axky - hy)^ . 

The valence bands in Eq. (HJi are described by the isotropic 
(72 - Ji) Kohn-Luttinger model, ^^ 



n^ 5 , 

Hv = -iz — [(ri + :^yi)k 



•2r2(k-S)2], 



(8) 



where nig is the electron mass, S represents the 3/2 spin op- 
erator matrices, y\ and 72 are two model parameters, and 
k - (kx, ky, k.) is the wave vector. For simplicity, we approx- 
imate {k,} - and {k^} ^ (^)^ in Hy when considering the 
quantization along the z direction of the quantum well, where 
d is the thickness of the quantum well. The eigenenergies and 
states of the valence bands are given by 



nhi 



^'72, 



ei/2(k) = e4/3(k) = -A - -^{k' + {kp) ± -^Q, (9) 



|l,k) 



13, k) = 



(Q] 




r ^ 




K 


, |2,k) = 


Q 



\ J 




I -K) 



( -K* \ 


Q 




|4,k) = 



f 

K* 



Q 



(10) 



where A is the energy gap, and Q - C^ (-k^ + 2{k^) + Q), 
K = V3Ci,feJ, Q. = 2Jk^ + <fe2>2 - kHkf), and C,, = 

3/^* + [k^ - 2{kl) - Qj" 

The electric-radiation interaction that couples the conduc- 
tion and valence bands is given by,^- 



V{t) = — A ■ p, 

llle 



(11) 



where -e is the electron charge and p the electron momen- 
tum operator. Under the electric dipole approximation, the 
vector potential A is related to the electric field of the single- 
frequency light by 



A(0 = — [E(w)e" 



■ E*(w)e''^'] 



(12) 



where we consider a linearly polarized single-color light inci- 
dent normally to the x-y plane, with the frequency w and the 
polarization electric vector 



E = £'o(cos6',sin6',0). 



(13) 



2m. 



2m 



where £0 is the amplitude and is the polarization angle with 
respect to the x axis. With the help of the Poynting vector, the 
electric component Eq can be evaluated from the energy flux 
of the laser by ^I - (1/2) ^JeqJJiqE^, where ^ is the absorp- 
tion efficiency and / is the light intensity, eo and //o are the 
dielectric constant and magnetic permeability in the vacuum, 
respectively. We always assume w > A, therefore the domi- 
nant optical absorption mechanism is the direct interband tran- 
sition. "Direct" means that the wave vector k of the electron 
keeps unchanged in the transition. 

We neglect the Zeeman effect in the valence bands, be- 
cause the contribution to currents from holes is expected to 
be much smaller than electrons considering the short charge 
and spin lifetimes of holes in n-type quantum wells. Besides 
the lowest subbands |-i-, k) and |-, k), another pairs of conduc- 
tion subbands |-i-', k) and |-', k) above them are also consid- 
ered in the numerical calculation. Approximated as the en- 
ergy levels of an infinitely depth potential well, |±') are about 
3{fP'/2m*)(7r/d)^ above |±). The differences in the effective 
mass and spin-orbit couplings between |+) and |+') are ne- 
glected for simplicity. Due to the parity in the z direction, 
conduction bands |+,k) and |-,k) couple to valence bands 
|1, k) and |4, k), while \+', k) and |-', k) to |2, k) and |3, k). For 
simplicity, we also neglect the diamagnetic contribution^^— 
which should give quahtatively similar results to the Zeeman 
effect. 



IV. ANISOTROPY OF PHOTOEXCITED CARRIER 
DENSITY IN k SPACE 

The photoexcited carrier density and all the physical quan- 
tities can be found within the density matrix formalism. In 
this work, we will consider only the steady-state nonequilib- 
rium photoexcited carrier density, which can be found by the 
approach similar to that of the second order nonlinear optical 
susceptibilitiesi^ We start with the Liouville equation of the 
density matrix. 



d,p„ 



-^[^'PL 



\Hnm Pnmj 



■y\p> 



(14) 



where n and m run over the states in Eqs. ^ and ( fTOl l. and 
y is a phenomenological damping parameter. p„^ is the equi- 
librium density matrix before the light excitation. Because we 
assume a «-type quantum well, the valence bands are fully 
occupied p', ^, = (5,,.,.', where v, v' e {1,2,3,4). While the ini- 
tial equilibrium density matrix of the conduction bands are 
described by the Fermi function. 



o(0) 



p:^;;:(k) = /[e,(k)]5,,,, = 



gfe(k)-£p]/<:B7- + I ■ 



(15) 



where c, c' € {+, -, +\ -'), E^ is the Fermi energy, k^ is the 
Boltzmann constant, and T is the temperature. By treating 
|c, k) and |v, k) as unperturbed part and V(t) as perturbation, 
the perturbation equations up to the second order are given by 



d,f^ - 


ICOnmPnm ' 






drP^l -- 


-- -i0J„,nP^;^l- 




-yf^:2r 


d,p2 -- 


- -il^nmPnm ~ 


T\yrp''' 

fl L Jnm 


-yt^. 



(16) 



where w„„,(k) - [e„(k) - e„,(k)]/^. p„,„ are functions of k 
because the light is momentum-free under the electric dipole 
approximation. After a straightforward derivation, the leading 
order of the light induced steady-state density matrix for the 
conduction bands is found to be of the second order: 






^[v„(k)-E(w)][v..,-(k) ■£*(«)] 



X [(1 - fc)6{u - Wcv) + (1 - /c')<J(w - Wf',,)] , 

(17) 

where the steady state is approximated by introducing the mo- 
mentum relaxation time t^,. 



p^^'(k)*T,5,p^^'(k,0 



(2), 



(18) 



Tg - Pmfn* je can be estimated from the mobility /i,„ and the ef- 
fective mass m* . Equation ( fTTl i recovers the result by employ- 
ing the semiconductor optical Bloch equationsi^ and Fermi's 
golden rules. ^" In the following, we will suppress the super- 
script of p*^^^ for simplicity. By substituting the velocity \cv in 
Eq. ( fTTl l by the position r^v according to 

<c,k|v;|v,k) = <c,k|4[r,-,//o]|v,k) = /w„,<c,k|r,|v,k), (19) 
in 



Eq. ( [TtT i can also be expressed as 

Pc,f'(k) = —- — > [w,.-„(l - fc)8{u - w„,) 
n'-uj ^-^ 

+w„,(l - /,0<J(w - «,-„)] [r„,(k) ■ E(w)][r„„(k) ■ E*(w)]. 

(20) 

By using the eigenstates in Eqs. ^ and ( fTOl i and considering 
the spatial wavefunctions of the eigenstates, the elements of 
the transition matrix Vcv = {xcv,yc\) — (Xvc^y'vc) ^re found as 



Xvc = {V, k|i|c, k) = Qc 



vf(%-0 (Q- 



V3 



)U* 



^(7^+ -J) -iK^^)W 



(21) 



and 



y^c = <v, k\y\c, k) = ioc 






-(Q + ^W 



(22) 



where we have defined the efl'ective dipole length 

flcv = <0,0|x|l,-l) (23) 

with |0, 0) and |1,-1) the spherical harmonic functions Fo.o 
and Fi-i, respectively. Oc, has the dimension of length. U, 
K, and Q in Eqs. ( |2T] ) and ( l22l i have been defined in Eqs. (|7| 
and ( fTOb . Note that the optical selection rules owing to the 
s- and /9-wave natures of the conduction and valence bands 
have been incorporated in Eqs. (l2Tl i and (l22l i. We neglect the 
density matrix of the valence bands because the charge and 
spin lifetimes of holes are much shorter than those of electrons 
for n-type quantum wells. 

In this work, we will retain only the diagonal part of the 
density matrix because when ignoring the broadening of the 
light frequency co, the coherent contribution from the off- 
diagonal part of the density matrix can be neglected. '*' The 
diagonal terms of the density matrix have clear physical mean- 
ing as the photoexcited carrier density. By substituting the 
polarization electric field vector E = £'o(cos0, sin^) into Eq. 
( |20|) . the carrier density excited to conduction band \c, k) can 
be written as the summation from different valence bands 
|v,k) 



e(k) = 2^Pcv,k 



(24) 



where p^x represents the carrier density excited from valence 
band |v', k) to conduction band \c, k) as a function of the wave 
vector k. p^k can be divided into three terms according to 
their dependence on the polarization angle 6, 



cos 



pcvM = p"v,k+p™: 



%cos20 + pf^^^sm26, (25) 



where 

9 2 I 



(1 - fc)6{cO - COcv)(\Xcv\^ - lycvf'). 




, ^, —(1 - fc)6(aj - Wcv)Re(2x„,3'„c). 



(26) 



In the polar coordinates (k^, ky) = k{cos (p, sin (p), we can trans- 
form the delta function of a> into that of k. 



9 2 / 

P°.. k = ^^^71^ a/ -d - fc)G[k,Aip)mk - k„.(ip)] 
n^ V £() 

9 2 / 

xd^nf - \yj\ 

« V So 
xRe(2A;cy}'yc), 



where we defined 



G[kcviip)] = 



\dajcvik,tp)/dk\ ^^^^^(^) 



(27) 



(28) 



G[A:e,,(^c)] is of dimension of second-meter"'. kcv{<p) is the root 
of w = (jOcvik, <p) for a given (/3, i.e., the wave vectors on the 
constant energy contour (see Fig. |2]i. The delta function 6[k - 
kci,((p)] restrict the values of A: on the constant energy contours. 
In Fig. [3j we show the calculated zero-field p° , ^, pl°\, and 
P'cvk ^^ functions of ip, for the carriers excited from 1 1 , k) to 
|±, k). It shows that p°^ j^ is an isotropic function of if, while 

p™'*!^ and pj'"|^ depend on cos Itp and sin 2if, respectively, in- 
dicating the total carrier density excited by the linear light 
is anisotropic in k space, p" , |^ are always positive and over- 
whelm pI"!"^"^ in magnitude, so the total photoexcited carrier 
density is physically positive. In addition, as indicated by the 
"max"s and "min"s, the carriers excited from the same valence 
band have different carrier density on the spin-split conduction 
bands \+,k) [Figs. Hal)-(cl)] and |-,k) [Figs. [3a2)-(c2)]. 
This difference is due to the splitting of the constant energy 
contours for |-i-,k) and |-,k) by the spin-orbit coupling, as 
shown in Fig. |2ja). Fig. [3]shows only the photoexcited carrier 
density from 1 1 , k) to | ±, k). Other pairs of conduction and va- 
lence bands also have similar anisotropic photoexcited carrier 
density in k space. 



Summarizing the above, the carrier density excited via di- 
rect interband transitions from valence band |v, k) to conduc- 
tion band \c, k) by the normal incidence of a linearly polarized 
light can be expressed as 

Pe.-.k = P° . ,„ + P™1 COS 2ip COS 261 + pf",„ sin 2ip sin 26, (29) 
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FIG. 3: At zero magnetic field (6,- = 0, By = 0), the calculated 
photoexcited carrier density prv.k in Eq. ( 127! ) as functions of the wave 
vector angle (f. Dark (light) represents positive (negative) values, 
"max" ("min") indicates the maximum (minimum) values. The p°^,^ 
term is always positive and overwhelms p^°^ • p'J'f'j;*'*'" are in units 
of (2;r^/T,e-a^,, -y'so/po/^^) x (second/meter), and have a dimension 
of meter"'. Parameters: T =77 K, Ep = 0.01 eV, co = 0.8 eV, 
m' = 0.04m„ yi = 11.97, 72 = 4.36, A = 0.764 eV, rf = 40 nm, 
a = 4.31 meV-nm, yS = 0. Numerical scheme is given in Appendix 

E 



where pd™'*'"" in general are functions of ip. In a weak in- 
plane magnetic field (say, less than 1 tesla), PcjIJ™''*" slightly 
depend on ip. In the absence of the magnetic field and the 
anisotropy of the spin-orbit couplings, pd™'*'^'" become inde- 
pendent on if. In this case, Eq. ( |29] l can be written as 



Pcr.k - Pcv 



+ p™' COS 2ip cos 26 + p™ sin 2ip sin 26, 



COS 



(30) 



where P(,,f°^'™ are constants of dimension meter" ^ 

The anisotropy of the photoexcited carrier density is the 
core of this paper. It can naturally account for the field and po- 
larization dependence of the magnetoelectric photocurrentS^ 
and the pure spin photocurrent3^ 



MAGNETIC-FIELD INDUCED ELECTRIC 
PHOTOCURRENT 
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A. Origin of co, c,, and Cy 

The electric photocurrent density along p(e {x, y]) axis can 
be found by summing the velocities of the photoexcited carri- 



ers 



■J]P"-'^'^c., 



Jit - ^ / jHcv,k^„,]^- 

',v,k 



(31) 
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FIG. 4: The product table of the photoexcited carrier density and velocity as functions of wave vector angle (p. Only the carriers excited from 
|l,k) to |+,k) are shown, [(a)-(c)] The zero-field photoexcited carrier density p+i'J;''^'". [(i)-(iii)] The zero-field velocity (i), and its variation 
with Bj- (ii) and By (iii), respectively. "(0)" stands for the values without the magnetic field, and "A...(i?v)" for the variation upon applying the 
field By. Dark (light) represents that the value is positive (negative), "max" ("min") indicates the maximum (minimum) values. The boxes 
mark the origins of the three terms in Eq. l|T). The parameters are given in Sec. IV Bl yS = 1 meV-nm is assumed. Numerical scheme is given 
in Appendix IBI 



From the eigenenergies Eq. (|6]l, the velocity along the x di- 
rection for conduction bands |±, k) can be found as 



1 de± ti , a 



(ak_, - hy) 



^ dk, m* ■ n ^(aky + h,,)^ + (ak, - h,)^ 



(32) 



Above and hereafter, we will ignore p in the analytical ex- 
pressions, and take it into account only in the numerical cal- 
culation. By rewriting k^ - k cos ip and ky - ksvcup in polar 
coordinates {k, ip) and expanding the velocities up to the linear 
order in h^^ and hy, we have 



sm2(p 



^'Ik ^ (— ^±t)cos<^+ -— - 
m* n Ink 



hx + 



sin^i^ 
hk 



hy. (33) 



i^j, 1^ in the current density formula ( |3T] ) is related to Eq. (l33T l 
by restricting k on the constant energy contours k^y. 

Rewriting the summation in Eq. (|3TT i into an integral in 
polar coordinates, putting into Eqs. ( |29] l and ( |33] |, and per- 



forming the integral over k, we obtain an integral over tp 



Jx = - 



p2;r 

V I d(pkcy,((p) 



i^nf^j. 



x[P° ..^ + Pc°; cos lip cos 29 + p™^ sin lip sin 29] 

,,h a sin 2^ sin^ ip 

x[( — kcv(ip) + c—) cos ip - c hx - c «v], 

TO* h 2hkcv(<p) hkcv{ip) ' 



(34) 



where ip s arctan(fc,7A:j), k = k\+ k^,, and kci,(<p) are the roots 

of w = (jJcviK ip)- In general, kcv{ip) and Pci™**'"" are functions 
of 1^ when B ^ andyS ^ 0. At this moment, if we ignore their 
ip dependence, the above integral can be readily performed. 



jABx,B,)^h,-^{Y'-pl) {"dip 
' (27!-)2 ^h Jo 



sin^^ 
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+''■'■ TTTI^y.^^""^^*" 20 I # sin 2,^ sin 2^, 
(2:t)2 ^ 2^! Jo 

+(^2^^-^''^^^"20. (35) 



The resulting three nonzero terms immediately recovers the 
form of the experimental fitting formula ([T]l. In this way, a 
clear relation between the anisotropy of the photoexcited car- 
rier density and the magnetoelectric photocurrent density is 
established. 

Above we ignore the ip dependence of p^!;™''''*'" and k^y re- 
sulting from finite B and ji. Their influences can be taken into 
account numerically. In Fig. |4]we present a product table of 
calculated v'j ^ andp^'J^''™ as functions of i^, for finite B and 
y6 (refer to Appendix iBl for the numerical scheme). The three 
nonzero contributions in Eq. ( [35l l are marked by the boxed 
sub-figures (iic), (iiia), and (iiib). Except the three nonzero 
terms, note that all the other sub-figures in Fig. |4]always have 
equal weight of positive and negative lobes. In other words, 
they yield zero when integrating over ip, and give no contri- 
bution to the electric current. All the pairs of conduction and 
valence bands have the same behavior, and all of them add up 
(actually cancel with each other) to give the total net electric 
current. 

In the analytic result ( |35] l, we have neglected the variation 



O,cos,sin 
cv'.k 

charge current j^ induced by By can be expanded as 



of p^;, 1^ ' when applying the magnetic field. In general, the 



}^,{B,) - -e2bo„,k(0)v^^(0) +p„,k(0)A4(B,) 



h 



:,i',k 



+Ap,,, k(B,)/^(0) + Ap„, k(Bv)A4(Bv)] + 0{Bl\ (36) 

where "(0)" stands for "at zero magnetic field", and "A. ..(By)" 
for the variation when applying the magnetic field By along 
the V direction. We already illustrated that the first term van- 
ishes, and the second term is consistent with the experiment. 
Besides, one can expect that the third term, i.e., the mag- 
netic field-induced variation of the photoexcited carrier den- 
sity, gives a contribution of the same order as the second term 
[see Fig. 13, and the last term is ignorably small. Therefore, 
we approximate the photocurrent density induced by the mag- 
netic field by 
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(Bv) - (/7-^^cp>,-('7T^X'^'°-^''>'^°'20 



+{n^Y^cp'^)K^m2e, (31) 
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where an extra "2t]" has been multiplied to account for the 
summation of the second and third terms of Eq. ( |36] |. 
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FIG. 5: The product table of the field-induced variation of the pho- 
toexcited carrier density and the zero-field velocity as functions of (f. 
The legend is the same as in Fig. |4] The boxed three terms, when 
added to the boxed terms in Fig. [4] gives the "2rj" factor in front of 
Eq. ( I37l (. All parameters are the same as those in Fig. |4] Numerical 
scheme is given in Appendix IbI 



B. Comparison with experiment 

Above we show that the photoexcited carrier density in Eq. 
d29] l explains the origin of cq, c,, and c^ terms in the fitting for- 
mula ([T) of the magnetoelectric photocurrent. Now we com- 
pute co/x/y, and see how close the minimal model can be when 
comparing with the experiment. 

We choose a set of parameters close to the experiment. 
From the experiment, we have that the temperature T =77 K, 
the Fermi level Ep - 0.01 eV, the light frequency a> = 0.8 
eV, the band gap A ^ 0.764 eV, the quantum well thick- 
ness c/ = 40 nm, and the Rashba spin-orbit coupling con- 
stant a - 4.31 meV-nm.— The Luttinger model parame- 
ters yi = 11.97 and 72 - 4.36 are adopted from those for 
Gao.47ln() ssAs,-*"* which has the similar components as those 
in the experiment (Gai_ JnjAs, x - 0.53 ~ 0.59 by graded 
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doping). The momentum relaxation time can be estimated 
by Te - ^^^-^ ~ 2 ps, with the effective mass m* ^ 0.04 
and the mobihty /i,„ ^ 84000 cm2/(Vs) at 77 K (about 7 
times larger than that at room temperature^^). The effec- 
tive dipole length Acv ~ 6. 7 A, is approximated by that for 
GaAs.'" Considering a reflectance of 0.3 and the absorption 
coefficients 9 X lO^cm"', we obtained the absorption efficiency 
^ = (1 - 0.3)[1 - exp(-9 X lO^^cm-i x 40nm)] k 2.5%. With 
the light power P - \5 mW and the light spot radius 5|/m, 
the light intensity is found as / x; 1.91 x 10'^ W/m^. The ex- 
periment observed that |c.i| + Ic,!; this may be due to a finite 
Dresselhaus spin-orbit coefficient yS, so that the spin-orbit cou- 
plings along the x and y directions are different, i.e.. Ax + /!,, 
since Xx - a - ji and /l, - a + ji. The value of/? is unknown 
to the experiment. In this following calculations, /? will be 
chosen as a variable parameter 
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FIG. 6: (Color online) Calculated j^ for different/?, (a) jxiG) at By = 
1 T (b) jAB,) at = 120°. (c) jAd) at B, = 1 T. (d) /,(S,) at 
6 = 120°. Parameters are given in Sec. lVBl 

In Fig. |6] we present the numerically calculated photocur- 
rent densities jx as functions of Bi, By and 6. We also compare 
the results for positive, zero, and negative/?. Qualitatively, the 
calculated results capture the main features of the experiment. 
Fig. |6jb) and (d) show that, the current densities are linearly 
scaled with the magnetic fields. At zero magnetic field, there 
is no current. Fig. |6](a) and (c) show that, the parallel field 
leads to only sin 20 dependence while the perpendicular field 
induces only the constant and cos 20 dependences. In sum- 



mary, the calculated current density can be formulated in the 
same form as Eq. ([T]i. 

Now we make some quantitative comparisons. In Fig. |6] 
the current densities are in units of ^ a; 0.09 A/m [see Eq. 
(IB2l i1. Fig. |6] shows that the dominant term is jx{By), consis- 
tent with the experiment, /v(fiv) can be as large as 8 x lO^^f 
when By = 1 T, i.e., about 0.72 x 10"^ A/m, comparable with 
the experimental estimate ~ 2 x 10"^ A/m.^^ The calculated 
current is smaller, probably because only hmited bands are in- 
cluded in the calculation. In tablelH we compare the calculated 
Co, Cy, and Cj with the experiment. The calculated c,, and Cj are 

TABLE I: Comparison between the theory and experiment for the 
parameters cq/I, Cx/I, and c,.// in formula (|T](. / is the light intensity. 
Co/x/y/I are in units of lO^'^'x (A/m)/(T-W/m^). Parameters are given 
in Sec.lVBl 



Co// 


Cjl Cx/I 


\cJCx\ 


Experiments^ 10.2 


-1.74 -0.94 


1.85 


Theory (/J/q- = -1/3) 2.2 
Theory 06 = 0) 1.9 
Theory (j6/q-= 1/3) 2.2 


1.6 -0.8 
1.2 -1.2 
0.8 -1.6 
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1 

0.5 



comparable with cq, while in the experiment cq is almost an 
order of magnitude larger than c,, and Cy . In Sec. |Vl we have 
shown that c,, and c , terms come from the photoexcited car- 
rier p™^'*'", which originate from the quantum interference be- 
tween two circular components of a linearly polarized light.— 
In contrast to the ideal situation assumed in the theory, the in- 
terference effect may be suppressed in the experiment, then Cy 
and Cx are reduced and cq is enhanced by a non-interference 
contribution. Besides, the sign of the calculated Cy is opposite 
to that of the experiment. This may be attributed to the dif- 
ference in band symmetries between the experiment and the 
model. Considering the simplicity of the minimal model, we 
expect that more sophisticated models may cover more reli- 
able details, but the minimal model is enough to offer a rea- 
sonable physical picture for the experiment. 



VI. ZERO-FIELD PURE SPIN PHOTOCURRENT 

A. Polarization dependence of zero-field Pure spin 
photocurrent 



r 



Both the symmetry argument of the C2v group'^ and theo- 
retical calculations "'■^"■^' have pointed out the generation of 
pure spin currents by the linearly polarized or unpolarized 
lights, and they can be expressed as functions of 6: 

jl = /o -H /, COS 26, il = h sin 26, (38) 

where /o,i,2 are constant coefficients. In this section, we will 



show how the anisotropy of the photoexcited carrier density is 
related to the ^-dependence of the pure spin current. 



The spin photocurrent density flowing along the ju (g {x, y}) 
direction with spin polarized along the v (e {x,y]) direction 
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FIG. 7: The product table of the photoexcited carrier density and spin velocities as functions of wave vector angle ip at zero magnetic field. Only 
the carriers excited from 1 1, k> to |+, k> are shown, [(a)-(c)] The photoexcited carrier density p+i_k- [(i)-(ii)] The spin velocities flowing along x 
axis, while with the spin pointing along x [(i)] and y [(ii)] directions, respectively. Dark (light) represents that the value is positive (negative), 
"max" ("min") indicates the maximum (minimum) values. The boxes mark the origins of the three terms in Eq. l |38t . The parameters are given 
in Sec. lVBl yS = 1 meV-nm is assumed. Numerical scheme is given in Appendix IbI 



can be found by 



J,- 



c.vM. 
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(39) 



where we also retain only the diagonal density matrix, and 
s^J^,^^ is the spin velocity for the carriers excited from valence 
band |v, k) to conduction band \c, k). The spin velocity opera- 
tor is defined as 



--5-4t' 



(40) 



where a-y is the Pauli matrix. The zero magnetic field expec- 
tation value of the spin velocity for conduction band \c, k), 
defined as /^^ = {+, k|.?'^''|±, k), can be found in polar coordi- 
nates as, 
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/^^,^^ in the spin current density formula ( [39] ) is related to Eq. 
(1411 1 by restricting k on the constant energy contours between 
valence band |v, k) and conduction band \c, k). 

Similar to the current density, the spin current density can 
also be rewritten into an angle integral. 
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Ik 



dipkcy{ip) f^likcv, tp) 



X [pn;^ + pIZ '^os 2^ COS 20 + p^^^ sin 2^ sin 2^] . (42) 

If we ignore the ip dependence of pc,'.™*''^" and kcvi^p), the above 
integral yields 

m) = ^lsin20j]c-^^?„p-, 
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(43) 



which give the 6 dependence in Eq. ( [38] l. In the presence of a 
finite y6, the result is not affected qualitatively, as shown in Fig. 
|7] where the nonzero contributions to Eq. (l38l l are marked by 
the boxes. 

B. Quick estimate of zero-field spin photocurrent from 
magnetoelectrlc photocurrent 

Because both the magnetoelectric photocurrents and the 
zero-field pure spin photocurrents originate from the same 
photoexcited carrier density, this allows us to find a relation 
between them. With the help of Eqs. ([37j and ( |43] l, the ra- 
tio of the zero-field longitudinal pure spin photocurrent to the 
electric photocurrent induced by the parallel magnetic field 
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FIG. 8: (Color online) Numerical comparisons between the magne- 
toelectric photocurrents and zero-field pure spin photocurrents. (a) 
jABy) X 250 at By = 1 T and ;j(0). (b) jAB,) x 250 at B, = 1 T and 
7^(0). The parameters are given in Sec. IV Bl For a direct compari- 
son, we neglect the -e and fi/2 in front of the current density and spin 
current density formulas. The above values of currents are in units of 
(2n^lTge-^aly)/{h^ ■ nm), and have a dimension of second" 'meter"'. 



turns out to be 
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(44) 



and the ratio of the transverse spin photocurrent to the perpen- 
dicular magnetic field induced electric photocurrent is given 
by 



fM 



jxiBy) 
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Note that 



These relations have clear physical meaning. 

v-i O,cos,sin O,cos,sin O,cos,sin nru:^ a:4¥^^ ^^ ^4^ +u^ 

2jf=+ cpev- — P+v ~ P-v ■ '^"^^ dirrerence or the 

photoexcited carrier density between the + and - conduction 
bands, as already shown in Fig. |3] is due to the spin-orbit 
coupling. Therefore, the denominators in Eqs. ( l44l i and (l45T l 
mean that both the magnetic field and the spin-orbit coupling 
are necessary ingredients of the magnetoelectric photocurrent, 
while the numerators indicate that the pure spin currents are 
proportional to the spin-orbit coupling and the kinetic energy 
of the photoexcited carriers (if we can view literally h^k^.^,/2m* 
as kinetic energy). 

At this moment, we make a bold approximation by cancel- 
ing the effect of the spin-orbit coupling from both the denom- 
inators and numerators in Eqs. (l44l i and (|45] l, and literally 
say that the ratio of the zero-field pure spin photocurrent to 
the magnetoelectric photocurrent is about "kinetic energy over 



Zeeman energy". This relation, though rather coarse, can help 
us to make a quick order-of-magnitude estimate of the unde- 
tectable pure spin photocurrent from the measured magneto- 
electric photocuiTent.^^ The average kinetic energy of the pho- 
toexcited carriers is higher than the equilibrium Fermi energy 
measured from the bottom of the conduction bands, and thus 
is more than Ep - 0.01 eV for our numerical calculations. The 
Zeeman energy induced by 1 tesla of magnetic field is about 
10 "* eV for the Lande ^-factor of g^ - —A^— Therefore, 
the rough estimate implies that the spin photocurrent is about 
two orders larger than the magnetoelectric photocurrent at 1 
tesla. To test the reliability of this quick estimate, we numeri- 
cally compare the magnetoelectric photocurrent and zero-field 
pure spin photocurrent in Fig. [8] The spin photocurrents are 
about 250 times larger than the magnetoelectric photocurrent 
at 1 tesla of magnetic field. Despite its roughness, this quick 
estimate gives a reasonable result. We expect that this quick 
estimate can serve as a reference for more sophisticated non- 
destructive approaches, such as the Faraday rotation^ or the 
second-order nonlinear optical effectSi^iiS 



VII. CONCLUSIONS 

In this work, we present a theoretical description of the re- 
cent experiments^ on the optical injection of spin-polarized 
carriers via direct interband optical excitations into a semi- 
conductor quantum well under the normal incidence of lin- 
early polarized or unpolarized lights. In that experiment, the 
injection produces pure spin photocurrents accompanying no 
electric current at zero magnetic field due to spin-orbit cou- 
pling that respects time-reversal symmetry. An in-plane mag- 
netic field can break time-reversal symmetry, and extracts a 
magnetoelectric photocurrent. The magnetoelectric photocur- 
rent is characterized by its dependence on the magnetic fields 
and the polarization of the light. 

With the help of the density matrix formalism, we calcu- 
late the photoexcited caiTier density, cuiTent density, and spin 
current density. The photoexcited carrier density in k space 
shows an anisotropic dependence on both the wave vector an- 
gle (/? and the polarization angle 6 of the linearly polarized 
light as given in Eq. ( l30b . Since the velocities of carriers 
can also be expressed as functions of ip [see Eq. ([33])1, we 
can show that the current density can be simplified as an an- 
gle integral over if of the product of the density and velocity 
of photoexcited carriers. The angle integral then produces all 
the magnetic field and polarization angle dependences of the 
magnetoelectric photocurrent as reported in Ref. |25l in partic- 
ular, the magnetoelectric photocurrent induced by the parallel 
magnetic field. 

We show that the present simplified model with Rashba and 
Luttinger Hamiltonians is able to reproduce the cuiTent for- 
mula for the magnetoelectric photocuiTent. However, discrep- 
ancies still exist between theory and experiment for the rela- 
tive magnitudes and signs among the parameters in the cuiTent 
formula. Since the simplified model is a natural choice con- 
sidering the symmetry of the quantum well that was investi- 
gated in the experiment, the discrepancies indicate that further 
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investigations, with more detailed considerations on the band 
structure of the sample, are needed to identify the origin of 
the discrepancies. For example, more energy bands may be 
required in the calculation, the graded doping in the quantum 
well may need a self-consistent calculation of the potential 
and carrier density. 

We further show that the origin of the previously predicted 
pure spin photocurrents^"-^' can be well illustrated from the 
same photoexcited carrier density. We propose that the ra- 
tio of the zero-field pure spin photocurrent to the magneto- 
electric photocurrent can be approximated as "kinetic energy 
over Zeeman energy". With this relation, the underlying pure 
spin photocurrent can be quickly estimated from the observed 
magnetoelectric photocurrent, and provides a reference for 
other approaches to nondestructive measurement of the pure 
spin photocurrent. 
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Appendix A: C2v group and anisotropic Rashba model 

It is well known that the heterostructures of inversion- 
asymmetric zinc blende materials grown along the [001] di- 
rection have the C2v point group symmetry. The basic build- 
ing block of these structures is shown in Fig. |9ja). It has four 
symmetry operations. When the x and y axes are defined as 
[110] and [110] crystallographic directions, respectively, the 
xz and yz planes coincide with the mirror reflection planes of 
the C2v group. As the basis functions, the polar vectors (such 
as velocity, current, electric field) along the x axis and the ax- 
ial vectors (such as spin and magnetic field) along the y axis 
transform according to the irreducible representation B i of the 
C2v group, and the polar vectors along the y and the axial vec- 
tors along the x directions according to the irreducible repre- 
sentation B2.^** The physical picture of "vectors transforming 
according to irreducible representations" is schematically il- 
lustrated in Fig. |9jb). 

The basic observation to the experiment data indicates that 
the current is linearly proportional to the magnetic field, and 
the 26 function dependence usually implies a second-order 
nonlinear optics. Phenomenologically, the current density j„ 
can be generally written as^^ 
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FIG. 9: (Color online) (a) The building block of an inversion- 
asymmetric zinc blende structure grown along [001] crystallographic 
direction. When x and y axes are defined along [lIO] and [1 10] crys- 
tallographic directions, the xz and yz planes coincide with the mirror 
reflection planes of the C2v point group, (b) How polar vectors /„ 
jy, E^, Ey and axial vector B^, By, cr^., CTy transform under the four 
symmetry operations of the C2v point group. E: identical; C2: two- 
fold rotation about z axis; cr,, ((t(,): mirror reflection with respect to 
xz (yz) plane. Whether the vector changes sign under the symmetry 
operation is indicated on the lower right corner of each panel by "1" 
or "-1", which are actually the characters of the representations Bi 
and B2 of the C2V group— 



electric field vector The nonzero terms of Eq. (lAll i require 
that the vectors on both sides transform in the same manner for 
all the symmetry operations of the €2^ group. Two examples 
are illustrated in Table |ll] In the language of the irreducible 
representation of group theory, the table |II]can be written as 



Bi ^ B2®Bi ®Bu 
Bi - B2®Bi ® B2. 



(A2) 



Similarly, all the nonzero terms can be found and summarized 
as Eq. ^. 

In addition, both cr^ and ky transform according to B 1 , while 
B 1 (g) B 1 yields the identity representation of the C2w group, 
so cr^ky is an invariant for structures with the C2v symmetry. 
Similarly, CTykx is also an invariant. On the contrary, o-jcky and 
CTyky are not invariants. As a result, the spin-orbit coupling up 
to the linear order in k can be generally described by the form 



Ja 



= X^P^'BpE.Ep, 



(Al) 



^SOC = /lvO"j:^y - '^xO-yky. 



(A3) 



where a, /3, y, and 6 stand for Cartesian coordinates. B^ and 
Ey are the components of the magnetic field and polarization 



where Ay and Ay are spin-orbit coupling coefficients along dif- 
ferent directions. In the coordinate system where x' II [100] and 
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TABLE II: Two examples of how to determine whether an element of 
the pseudo tensor ;f"'''"' is nonzero. 1 and -1 correspond to "remain- 
ing unchanged" and "changing sign", respectively, upon applying the 
symmetry operations of the Ciy group to the vectors, (a) For;^^^^^^, 
which is zero because /,- and B^E^E^ are different for cr,, and a[. (b) 
Por;^''''^', which is nonzero because j^ and Bj^E^Ey are the same for 
all symmetry operations. 
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y 1 1 [010], the spin-orbit coupling in the conduction bands of a 
sample with the C2v symmetry can be generally written as^^i^i 



"soc 



a{a-^.ky, - a-yk,,,) +P(o-,.k,, - cryky,) (A4) 



referred to as Rashba {a) and DresseUiaus {fi) terms, respec- 
tively. While in this work the coordinate system is ji(;||[110] 
and ^[[[llO], the form of the spin-orbit coupling can be ob- 
tained by rotating the above H'^qq by 45° 



^soc = (a+/3)cr_,ky-(a-/3)o-yk,. 



(A5) 



The extra /3 means that although the C2v symmetry allows 
only cr^ky- or cr,,A;j-type spin-orbit coupling when x||[l 10] and 
yl I [ 1 1 0] , the Dresselhaus term when x' 1 1 [ 1 00] and y' 1 1 [0 1 0] can 
lead to an anisotropy of the spin-orbit coupling. Usually, /3 is 
smaller than a?^'^^ The anisotropy of the spin-orbit coupling 
may explain why Ci ^ c,. in the experiment.— Note that this 



result applies only for when x and y directions are referred to 
the [110] and [110] cry stallographic directions. ^*^ 



Appendix B: Numerical scheme 

For each pair of conduction (c) and valence (v) bands, we 
denote all the quantum states by the discrete values of the 
wave vector angle tp ^> (pi - iln/N, with / = 0, 1, ..., A^ - 1. 
At each (/),-, the energy conservation law w - ojcvik, (fi) is nu- 
merically solved. The root of k is the wave vectors kcvifd 
on the constant energy contours. The physical quantities such 
as v'^^, y. and P^J^^^'^'" then are calculated by giving them k = 

kcvi(Pi){cos(pi, simpi), and denoted as vJ"^,^((/3,) andp°;™''""((/;,). 
As an example, we show how to numerically calculate the 6- 
independent term of jx- According to Eqs. ( l34b and dZTJ i. it 
means that we calculate 



/■I 3 2 I — ; — N-l 



X {l-Mkcr(fd]] 



c,v /=0 

;vP + lyj' 



G[k,,(<fii)], 



K-AiPi) 



2n J 



where Akcv((pd = 4kcv{ipi) 



\Kvi<Pi) 



(Bl) 
v(^/-i)|, and 



2™ y ^l<^cv(fi) gives the circumference of the constant energy 
contour In the summation, we are using the units [G] ~ sec- 
ond/meter, [v"] ~ meter/second, [Akcy] ~ nm"', so the units 
of the current density in front of the dimensionless summation 
of Eq. dBlb are 



[;°] ~ \ " * 0.09 A/m. 



(B2) 
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